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Abstract 

We analyze generation of phonons in tunneling structures with two electron states coupled by 
electron-phonon interaction. The conditions of strong vibration excitations are determined and 
dependence of non-equilibrium phonon occupation numbers on the applied bias is found. For high 
vibration excitation levels self consistent theory for the tunneling transport is presented. 



Tunneling current induces generation of phonons (or vibrational quanta for a molecule) leading to 
effective "heating" of the phonon subsystem in any real system with electron-phonon interaction. In 
scanning tunneling microscopy experiments this effect may induce motion, dissociation or desorption of 
adsorbed molecules thus allows single molecule manipulation on a surface 0,0, 0- The main purpose 
of the present work is to reveal the conditions for strong phonon generation as well as for its suppression. 
We also investigate the influence of strong phonon generation on tunneling conductivity behavior. 
In our previous work we analyzed modifications of the tunneling current using a model in which 
O 'electron-phonon interaction leads to transitions between two electron levels. In the present paper we 
^ 'demonstrate that in this system tunneling current can induce strong phonon generation. The mechanism 
q .of this enhanced phonon heating is closely connected with the presence of at least two electron states 
in a molecule or a quantum dot, and in a single level model, widely discussed in literature (0,0, 0), 
this mechanism is absent. 

We consider a tunneling system, which is described by the Hamiltonian of the following type: 
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H = Hdot + H tun + H (1) 

The part H do t corresponds to a QD or a molecule with two localized states taking into account. Electron- 
phonon interaction leads to transitions between these two states. 

Hdot = £ i a t a i + g( a t a 2 + ajai)(b + b + ) + u b + b (2) 

i=l,2 

where e$ corresponds to discrete levels in quantum dot (or two electron states in molecule) and we adopt 
that e\ > 62 i — optical phonon frequency (or molecule vibrational mode)and g — is electron-phonon 
coupling constant. Tunneling transitions from the intermediate system are included in H tun 

H tU n= T Pti (c+ai + h.c.) + Y T k,i( c £ a i + h-c.) (3) 

p,i=l,2 k,i=l,2 

And free electron spectrum in left and right electrodes (k and p) includes the applied bias V: 

H o = J2( £ k ~ A*) c k c k + J2( £ p ~ V - e ^) c p c p ( 4 ) 



* e-mail: ars@lpi.ru 



1 



Operators Ck, c p correspond to electrons in the leads and eij - to electrons at the localized states of 
intermediate system with energy e,. 

By means of Keldysh diagram technique non-equilibrium phonon numbers can be found from Dyson 
equations for phonon Green function. To determine phonon occupation numbers we need to solve the 
Dyson equation for D < together with D R ^: 

D<(ft) = D<(ft) + D<(Q)U A (Q)D A {Q) + D R {Q)U < (Q)D A (Q) + D R (Q)U R (Q)D < (Q) 
D R (n) = D R (ty + D R (£l)Il R (n)D R (ty (5) 

where Dq is equilibrium phonon Green function: 



(ft + id) 2 - ujI 

D<(Q) = N (Q)(Do^)-D A m 

Nq is Bose distribution function. Polarization operators in the lowest order in electron-phonon interac- 
tion (first order in g 2 ) are easily determined: 

U A (Q) - U R (Q) = -Aig 2 [ ^ \lmG A {uj)ImG A (u - tt)(n 2 (uj) - n x (uj - ft)) 
+ ImG A (uj)ImG A (uj - ft)(m(u;) - n 2 (u - ft))" 

n<(ft) = -4ig 2 [ — \nAu)(n 2 {u - ft) - l)ImG A (u)ImG A (u - ft) 
J 2n L 

+ n 2 (oj)(n 1 (uj-n)-l)ImG A (uj)ImG A (uj-n)\ (6) 

All electron Green functions in the above expressions are calculated with full account for tunneling 
transitions. Thus electron non equilibrium filling numbers nx,n 2 are determined by the tunneling 
processes (neglecting electron-phonon interaction) as: n^uj) = {^n° p {u) +^n° k {oj)) /^i. Tunneling rates 
7 are determined as usually by the corresponding tunneling matrix elements I^ fo) i and densities of 
states v kp of the leads: 7* = T ki i/ k , 7? = T pi v p . Total width of electron levels due to the tunneling 
coupling is denoted by 7« = 7f + 7 4 fc . 

From Eqs.(0) one can easily derive D < : 

D< : = n^iv* '- d a ) (7) 

Equation © together with the relation D<(Q) = N(Q)(D R (Q) - D A (Q)) leads to the following 
nonequlibrium phonon numbers: 

n< 

jV(ft = — _ (8) 

Polarization operator Il < can be easily divided into two parts: 

n<(ft) = 2iImU A N (n) + iP<(ft) (9) 

Which allows explicitly separate the equilibrium distribution function and nonequilibrium occupation 
numbers: 

iV(ft) = iVo(ft) + AiV(ft) (10) 

Where 

AiV(ft) = (11) 
v ; 2ImU A v ; 



Substituting in Eq.© for ni,n 2 their explicit expressions in terms of Fermi functions in the leads n p 
and n k we obtain: 

p<(n) = [ du (n° p (oj - n) - n° k (u; - fi)) 



7i72 



x 



And 



JmGf (u)ImGZ(u - n) {lH<^) - 7fryfn°(a;) + (7i7! " 7 2 Wo(^) 
+ JmGf ( W - n)/<( W ) (7MM - 7h^(c) - (7*7* - W) 

7i72 



'12) 



////I .I- 1 = — — / du [lmG A (co)ImG^(uj - tt) x (13) 

( 7l fc 7 2 (n°(o;) - n°(c - fl)) + 727?KM - n° p (u - SI)) + ( 7l & 7 f - lh P i)K(uj - n) - n°(« - fl))) 
+ImGf (w - Sl)ImG 2 A (iu) x 

(7 2 fe 7i(n fc H - ^ - n)) + 7i7?K(^) - n°(a; - Q)) - (-yf-yf - tItiOK^ - A) - n°(a; - n))} 

In the previous paper (jlj) it was shown that two different types of inelastic tunneling current 
behavior exists, depending on the ratio between four tunneling rates 7i,7i,7 2 7 2 - The sign of the 
combination (7^72 — 7 2 7i) determines the relative population of the two electron levels due to the 
tunneling current through the system, because the following relation holds: 



k P k P 

75 [72 ~ 7j7i 
(7172) 



(n° k (uj) - n° p (uj)) = m(u) - n 2 (uj) 



= 1 

8 2 = 2 
w=0.4 

g=0.04 



y = 0.08 



Y" =0.12 



y; = 0.14 



y 2 = 0.08 



Considering phonon generation processes cases of "normal" (ni(u) — n 2 {uo) < 0) and inverse (ni(u) — 
n 2 {uj) > 0) population for the two levels with e\ > e 2 drastically differ from each other. In the 
case of normal occupation phonon generation is rather weak. Typical dependence of phonon occu- 
pation numbers on bias voltage calculated by the help of Eq. fjllj) is depicted in Fig.l. Maximum 
value of nonequilibrium phonon numbers does not exceed several units. Increasing of temperature 
smears fine structure of this dependence and always decreases maximum values of phonon generation. 

On the contrary in the case of inverse occupation at some 
bias voltage the nonequilibrium generation infinitly in- 
creases. The divergence of nonequilibrium phonon fill- 
ing numbers occurs because at some bias voltage ImIl A 
given by Eq.([13|) passes through zero, changing its sign. 
In order to describe correctly the threshold of phonon 
generation it is necessary to take into account higher or- 
der terms which allows to take into account nonlinear in 
phonon occupation numbers effects. For "normal occu- 
pation" case ImIi A is always positive for any voltage and 
the Eq.([ll|) is sufficient to calculate the non-equilibrium 
phonon filling numbers, 
v In a tunneling system with inverse electron level 

occupation one should self-consistently take into ac- 
Figure 1: N p h m on versus applied bias for weak coun t modifications of electron Green functions due to 
generation regime. Solid line corresponds to e lectron-phonon interaction together with nonequilib- 
zero temperature, and dashed line - to T=0.2 r j um phonon numbers. The first corrections to the elec- 
tron Green functions and to the electron-phonon vertex of the order of g 2 D(Sl) (where D is dressed 
phonon Green function) is shown in Fig. 2. 
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Figure 2: Diagrams for n« 



Diagrams in Fig.2a describe the part of ImH^ A , which is connected with the corrections to electron 
Green functions. They can be analytically expressed by Eq.fjBJ), with g 2 order corrections either to 
ImGf, or to ImG A , or to occupation numbers n\ and n 2 . Let us denote this part by ImIl^ A 



(14) 



iImIi {l)A (9) = -2tg 2 [ ^ \lmG 2 1)A (u)ImG A (w - Q)(n 2 (uj) - n^u - SI)) 

y J Z7T L 

+ ImG A (u)ImG { l )A (uj - Q)(n 2 (uj) - n x {u - Q)) 

+ ImG^{u)ImG A {uj - Q)(n (1 \uj) - n (1 \uj - Q)) 

+ (1< — >2) 

Both the changes in electron spectral function and in occupation numbers are determined by self-energy 
parts Y,(lo): 



ig~ I [D R {uo')G<{uo - uo') + D > {u')G*{u - u/)]du/ 



S<(cj) = -ig 2 J D < {u')G<{u-u')du > 



(15) 



(for G 2 is replaced by G±.) 

From the Dyson equation corrections to electron occupation numbers are: 



n£\w) 
n 2 l \uj) 



|S<(^) + n 1 (a;)JmSf(cu) 
7i 

|S 2 < (cj) + n 2 (uj)Im^{uj) 



(16) 



72 



Self-energy parts S are determined by the Eqs.()15|). where phonon Green functions depend on non- 
equilibrium phonon numbers, which should be determined self-consistently. As we shall see below, in 
the weak coupling regime the width of the phonon lines remains narrow enough even in the presence 
of strong phonon generation. This allows us to use 5-function approximation for ImD R throughout all 
the calculations. For and n 2 we obtain: 



n {X \uo) 



n 2 1] (uj) = —ImG A x {uj - uj ) [n 2 (w)(ni(a; - uj ) - 1) + N(w )(ni(w - u ) - n 2 (uj))} - (u < — ► -u ) 

72 

It is important that some corrections appeared in ImH^ A are proportional to the nonequilibrium 
phonon numbers N(uq). These terms ensure nonlinear limiting of phonon generation. In the following 
we retain only these terms, because the rest part of ImH^ A , independent of N{uq), gives only very 
small shift of threshold voltage for strong phonon generation and can be omitted. 



7i 
9 2 



ImG A 2 (uj - ujq) [ni(uj)(n 2 (uj - cu ) - 1) + N(u )(n 2 (u - cu ) - ni(a>))] - (tu 



-Uo) 
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In the present problem vertex corrections to polarization operator shown in Fig. 2b are not so 
small as it is usually supposed for bulk electron phonon interaction, and also should be taken into 
account. Collecting the terms with phonon occupation numbers for diagrams in Fig. 2b. we get the 
second correction to ImIi A proportional to N which we denote by ImIi^ A : 

duo 



x 
2tt 



/mnW^(O) = -2g 2 N(u; ) |/ 

x Re [G A {uj)G a {uj - u )(Gf(uj -Q- uj )G*{uj - O) + G a {uj -Q- uj )G<(uj - Q)) 

+ (G A (u)G<(u - u ) + Gf(u)G* (u - uj ))G?(u -Q- u )G$(u - Q) 

+ G a {u)G a {uj + u )(Gf(uj - Q + u )G* (u - ft) + G a {uj -Q + uj )G<(uj - Q)) 

+ (G A (uo)G<(uo + ujq) + Gf{uj)G* (u + u ))G? (u -Q + lu )G*{lj - Q) 

+ (1"2)} (17) 

Integrating Eq.fjHJ) over Q near u Q we get the equation for self-consistent calculation of non equilib- 
rium phonon occupation number N p h(uJo) — I dQN(Q)ImD A (Q): 



2i(ImU A (uj ) + ImUW A (uJo)) v ' 

where ImU {1)A = ImU^ A + ImU^ A . There is no need to calculate corrections to the function Il < (a;o), 
because it is always positive and never becomes close to zero. So all the highest order corrections are 
inessential in this case. 

Finally we can rewrite Eq.(|18|) in the following form: 

1 0) 2ImU A (uJo)+g 2 N(uj )A(uj ) { ] 

Where P < and ImYi A are determined by the equations (jllj) and (J13)) . and all the essential second order 
corrections to ImIl A proportional to N(uq) are now rewritten as g 2 N {ujq) A{ujq) . The function A{ujq) 
depends on applied bias voltage through electron filling numbers in the leads. Cumbersome expression 
for A(ujq) is presented in the Appendix for completeness. 

Equation (fTU|) is a simple quadratic equation with the solution: 



g 2 A(uj ) L 



+ J(ImU A (u )) 2 + (N (u )2ImU A (u ) + P<(lu )) 9 2 A(uj ) 



(20) 



For normal electron level population ( g 2 A(cuo) ImU A (uo) in all bias range) Eq. (j20)l gives only 
small corrections to the first order Eq. (jllj) for nonequilibrium phonon numbers. 

For inverse population the appearance of correction g 2 A(tuo) is crucial, because ImH A (uJo) changes 
its sign and is equal to zero at some threshold value of applied voltage. It was found out that function 
A(ujo) also changes its sign (sometimes not once) being positive at large bias. Moreover the bias value 
at which A(ujq) goes from negative to positive values last time is very close to the threshold value for 
ImU A (u ). This is the reason why we should retain all the corrections of the second order to polarization 
operators in order to get the reasonable accuracy of calculations. It was checked by direct numerical 
investigations, that for some parameters of the contact only vertex corrections ensure the validity of 
Eq. (j20|) in all bias range. 



For large voltages, greater than the threshold value, ImIl A becomes large negative value. So we get 
from Eq.(|2U|)the following expression for phonon occupation numbers: 

For voltages beyond the threshold and if |ei — e 2 — uj \ > 71 + 72 the value of ImU A (u ) can be estimated 
as: 

|/mrr>o)l * -^-{iH - iH) ( 7 72 , 2 (22) 

7! 7 2 (e 1 -e 2 -uJo) 
For resonant case the denominator (e\ — e 2 — uj ) is replaced by (71 + 72) in this expression. The value 
of A in the range of large voltages can be estimated retaining only the following resonant term of the 
total expression: 

A{uo) = -8g 2 I p(ImGf(u;)) 2 (ImG A (co - u )) 2 (n 2 (u; - cu ) - m(u;)) + - (23) 
This gives for saturation regime: 

A M * JL -» + ^ (24) 
7i72 (£1 -£2 -0JqY 

Therefore maximum occupation numbers in saturation regime at high voltages are: 

N ^ (iW- ihV (25) 
g 

Since we consider the case of weak electron-phonon interaction, g <C 71, 72, phonon occupation numbers 
are large, which means strong phonon generation. 

Now let us return to the problem of self-consistency of the presented second order calculations. The 
broadening of electron levels due to electron-phonon interaction is determined by ImT> A fJ15j) . Our 
approximation remains valid until ImT, A is less, than the broadening of electron levels 7^2 due to the 
tunneling coupling. For large phonon occupation numbers the main term of ImY> A is: 

ImSf(w) = g 2 N(uj )(ImG A (uj - u ) + ImG A (u + u )) + ... (26) 

Thus N(u Q ) is limited by inequality: 

a 2 

maxImZf ~ — N{u ) < 71 (27) 

72 

From Eq. (|2fi|) the limit of validity of our approximation is: 

maxN < ^ (28) 
9 2 

This value of N max from Eq. (j23|) is just of the same order. So at large values of applied bias we, strictly 
speaking, work at the limit of applicability of suggested scheme but never go beyond this limit. 

The phonon line width is determined by the imaginary part of the polarization operator U A . Omit- 
ting second order corrections, maximum value of ImU. A ~ (? 2 /(7i + 72) (see ([22))) is much smaller than 
the phonon frequency uq. Second order corrections to ImIi A do not change substantially this value, 
even for large nonequilibrium phonon numbers. In accordance with Eqs. (126125(1 maximum possible 



changes of electron Green functions are of the order of themselves. Thus corrections to polarization 
operators can result only in some numerical coefficient of the order of unity. This means that even in 
the saturation regime with large nonequilibrium phonon numbers the width of the phonon line remains 
small. 

Tunneling current through two-level system can be expressed as jl],jH] : 

I(V) = 2 £ 7* / [ImGf (u)n° k (u) - iG<{u)]du (29) 

Corrections to the current are thus determined by the corrections to the electron Green functions G^ 2 
and Gf 2 . First order corrections corresponding to the first diagram in Fig. 3 with equilibrium phonon 
Green function were calculated in the previous paper j3] . If we take into account nonequilibrium changes 
of the phonon occupation numbers, we should consider four diagrams, depicted in Fig.3. 



Figure 3: Corrections to electron Green function 

It is worth to remark that just these types of the electron self-energy corrections are consistent with 
phonon polarization operators, depicted in Fig.2. The both sets of diagrams originate from the the 
same generating functional. 

This self-consistency means, that corresponding Ward identities are satisfied, and charge conserva- 
tion in the tunneling processes is automatically fulfilled. Note that all electron Green functions in all 
polarization operators and self-energy parts are calculated omitting the electron-phonon interaction. 
An attempt to use self-consistent Born approximation for the electron GF leads to violation of charge 
conservation in this problem and artificial symmetrization is needed to restore current continuity 

ma)- 

Nevertheless we can neglect the contribution from the three last diagrams in Fig.3. in the following 

cases: 

a) If there is no inverse electron population due to the tunneling current, then phonon numbers are 
small, and second order diagrams have small parameter g 2 /jij^. 

b) If inverse population appears, but phonon frequency is far from the resonance between electron 
levels: \e\ — e 2 — c^ol > 7i,72- Then second order diagrams have small parameter 7i/|ei — s 2 — ujq\. In 
this case phonon numbers are large N ^> 1, so the first correction with full phonon function strongly 
differs from equilibrium result. 

c) If inverse population appears and phonon frequency is almost in resonance with electron transition 
energy \e\ — e 2 — ujq\ > 71,72, then we can restrict ourself to the first term only near the threshold 
voltage, until phonon numbers are not too large. For high voltages in saturation regime in this case we 
can not use perturbation theory, because corrections to electron Green functions become of the order 
of zero order Green function. We should sum up not only second order diagrams depicted in Fig.3, but 
all the higher order terms as well. 

So we calculate the corrections to the tunneling current using only the first term in Fig.3 ( one-loop 
diagram) but with full nonequilibrium phonon green function, taking into account that in resonant case 
c) it makes sense not far beyond the threshold of strong generation. 

In addition to the tunneling current calculated in jl] there exists a correction to Eq. (}2*9"|) proportional 
to nonequilibrium phonon numbers: 



a r r k p 

AI = 2g 2 N(u ) [ ^Im{(G^ )A (iu)) 2 }(ReGt(uj-uo) + ReG^(uj + iu )) 
J 2,71 71 
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+ Re{(G^ A (uj)) 2 }(ImG$(uj - oo ) + ImG^oo + u ))] (n° k (uj) - n°Joj)) 



7*72 ~ 7*7l 

7i72 



{ni(u) — n2(uJ — uJo))ImGf(uj)ImG2(u — ooq) + (1 < — ► 2) 



(30) 



Since N(uo) depends on applied bias increasing rapidly at some threshold voltage, new peculiarities 
connected with dN/dV can appear in the tunneling conductivity (dl/dV). 

This peculiarity is more pronounced if u > {e\ — e 2 )- The peak in dl/dV at phonon generation 
threshold voltage (when dN/dV ^> 1) is clearly seen in Fig. 4a. 



10 



81 12 T k =0.11 7^=0.07 



e 2 = 1.5 

03=0.7 

g=0.04 



Y 2 =0.08 Y 2 =0.12 



— > 1 ' 1 > r^-^ r 

0,0 0,5 1,0 1,5 2,0 2,5 



(a) 



3,0 

V 




V 




Figure 4: N p h on on and dl/dV versus applied bias for strong generation regime. In the left panel (a) 
conductivity without electron-phonon interaction is shown by the dashed line. In the right panel (b) 
strong generation regime with non- monotonous behavior of N p h on on is shown. Correction to dl/dV 
(dotted line) is pronounced but can not be resolved in the total conductivity (solid line) 



For Uq < (si — 62) nonequilibrium phonon numbers can depend on applied bias non-monotonically 
(Fig.4b.). But in this case, in spite of great phonon excitation, no definite peculiarities connected with 
electron-phonon interaction can be observed in tunneling conductivity curves (see Fig.4b.). 

Approach to saturation regime for iV at large bias is clearly seen in Fig.4. The value of N, 
is in agreement with the estimation given by Eqs.(jz 



The obtained values for N„ 



max 

30 - 40 



for reasonable junction parameters can hardly be observed in real small systems, because such great 
overheating should lead to dissociation, desorption of molecules or destruction of quantum dots. In 
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order to deal with real objects at such great intensity of vibration excitation one should take into 
account phonon anharmonicity and other relaxation processes for phonons. 

In conclusion we point out that intensity of phonon generation can be tuned by changing the pa- 
rameters of the tunneling junction (which influence the tunneling rates). For inverse population of 
two-level electron system strong phonon generation takes place which can lead to drastic changes of 
the properties of nanostructures. Let us note that the problem of suppression of tunneling current 
induced phonon generation is very important for fabrication of semiconductor cascade lasers based on 
sequence of tunneling junctions [H]. Generation of optical radiation requires the inverse population of 
two electron states. But, as we have shown in the present work, in this case strong phonon generation 
inevitably appears and always competes with optical radiation generation. Using the results of the 
present paper one could analyze whether it is possible to achieve the threshold of optical generation 
before strong phonon generation begins changing the parameters of the tunneling system. 

This research was supported by RFBR grants and 04-02-19957, grant for the Leading Scientific 
School 4464.2006.2 and RAS Program "Strongly correlated electrons in metals, semiconductors and 
superconductors" . Support from the Samsung Corporation is also gratefully acknowledged. 
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1 APPENDIX 

Complete expression for the function A in Eq.(|19|) is the following: 



A(u ) = -4g' 



^2(/mGf (uj)) 2 (ImG A (cu - u )) 2 (n 2 (u - u ) - m(w)) 

Z7T 

1 



+ ImG A {uj)ImG 2 (uj - uj Q )ImG A (uj - 2u Q )(ImG A (u - uj ) )(n l (u - 2lu ) - ni{u)) 

72 

+ \lm{{G^{uo - uo Q )) 2 ){ReG A {uo) + ReG A {uo - 2u )} x 

x [ImGf(u>)(ni(u>) — n 2 (oj — u> )) + ImG A (tu — 2lu )(ti 2 (uj — lu ) — t%i{u — 2uj ))] 
+ ImGf(uj)ImG A (uj — ujq) ImG 2 (uj — uJo)ImG A (uj — 2iOo)(^n 2 (uj — lo ) — ni(uj) — rii(uj — 2oo )) + 

+ ImG 2 (to + tu )ImG 1 (co)(n 2 (uj — cu ) + n 2 (cu + tu ) — 2ni(u)) 

+ ReG A {uj)ImG 2 {uj — uq) ReG 2 (u — u )ImG A (u — 2uo)(n\(uj — 2uo ) — n 2 (uj — oo )) + 

+ ImG 2 {ijj + oo )ReGf(uj)(n 2 (uj + uj ) — n 2 {uj — oo )) + ReG 2 [uj + uj )ImG A (u)(rii(u) — n 2 {uj — lu )) 

+ ImG A (uj)ReG A (u - u ) ReG A (ui - u )ImG A (ui - 2u ){n 1 (u) - n^u - 2lo ))+ 

+ ImG 2 {uj — uJo)ReG A (u — 2uj )(rii(uj) — n 2 {uj — ujq)) + ImG 2 {uj + uJo)ReG A (u)(ni(uj) — n 2 {uj + ujq)) 

+ ImGf(uj)ImG 2 (uj - lu ) \ReG A (u - uj )ReG A (uj - 2oj )(n 1 (u) - n 2 (u - lo )) + 

+ ReG A (u + ujQ)ReG^{uj){ni{uj) - n 2 (cu - u )) 

+ ReG A (u)ReG 2 (uj — u ) ImG 2 {uj — uJo)ImG A {uj — 2uj )(n 2 (uj — u ) — ni(u> — 2uj )) + 

+ ImG 2 (u + uo)ImG 1 {uj)(n 2 (uo + ujq) — ni(uj)) 

+ (1< — >2) (31) 

Typical example of polarization operators behavior for the case of normal electron population is shown 
in Fig.5. 
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Figure 5: Dependence of polarization operators n < , ImIl A and A on applied bias for the case of normal 
electron population. (e 1 = 1.2, e 2 = 1.5, u = 0.4,0 = 0.04, 7l fc = 0.08,7* = 0.1,7? = 0.09, t| = 0.08) 

ImIl A is always positive and large enough, so corrections included into the function A are inessential. 
But for the inverse population we observe quite different behavior of ImIl A , depicted in Fig.6. If we 



in 




Figure 6: The same polarization operators for the inverse population. In the right panel enlarged region 
near the bias, where ImU A passes through zero, is shown. Blue line - complete function A, orange line - 
vertex corrections, green line -corrections only to electron Green functions (ex — 1, £2 = 2, u — 0.5, g = 
0.04, 7l fc = 0.14, 7 2 fc = 0.08, 7? = 0.06, 7! = 0.12) 
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